Quasi regular concentric waves in heterogeneous lattices of coupled oscillators 
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We study the pattern formation in a lattice of coupled phase oscillators with quenched disorder. In 
the synchronized regime concentric waves can arise, which are induced and increase in regularity by 
the disorder of the system. Maximal regularity is found at the edge of the synchronization regime. 
The emergence of the concentric waves is related to the symmetry breaking of the interaction 
function. An explanation of the numerically observed phenomena is given in a one-dimensional 
chain of coupled phase oscillators. Scaling properties, describing the target patterns are obtained. 



The study of coupled oscillators is one of the funda- 
mental problems in theoretical physics and has led to 
many insights into the mechanisms of spatio-temporal 
pattern formation in oscillatory media 0,12> > with ap- 
plications in a variety of systems such as arrays of Joseph- 
son junctions or oscillating chemical reactions Re- 
newed interest stems from its possible role in many bio- 
logical systems like cardiac tissue gg, neural systems |jj 
and ecological systems However, nearly all theoreti- 
cal studies have been carried out with idealized systems 
of identical oscillators, whereas not much is known about 
the dynamics and pattern formation in heterogenous os- 
cillatory media Here we show that spatial disorder 
leads to the emergence of quasi regular concentric waves. 

Target waves are one of the most prominent patterns 
in oscillatory media and are usually associated with the 
presence of local impurities or defects in the system []J, 
lit!. lllj . These pacemakers change the local oscillation 
frequency and are able to enslave all other oscillators in 
the medium, which finally results in regular ring waves .1, 
ll^ l . However, the assumption of a discrete set of localized 
pacemaker regions in an otherwise homogeneous medium 
is somewhat artificial. Especially biological systems are 
often under the constraint of large heterogeneity. In such 
a disordered system no point can be distinguished as a 
pace-maker and it is not clear whether such a system can 
sustain highly regular target patterns and where they 
should originate. 

In this paper we investigate the influence of quenched 
disorder on the pattern formation in a lattice of coupled 
phase oscillators. As we show the random nature of the 
medium itself plays a key role in the formation of the 
patterns. As the disorder in a rather homogeneous syn- 
chronized medium is increased we observe the formation 
of quasi regular target waves, which result from an intri- 
cate interplay between the heterogeneity and a symmetry 
breaking of the coupling function. 

We study a system of N coupled phase oscillators 0] 



I = 1, 



(1) 



one- or two-dimensional lattice. The natural frequencies 
u>i are fixed in time, uncorrelated and taken from a dis- 
tribution p{u>). A scaling of time and a transformation 
into a rotating reference frame can always be applied so 
that e = 1 and the ensemble mean frequency lJ is equal 
to zero. We refer to the variance a 2 = var(cji) of the 
random frequencies as the disorder of the medium. 

The effects of coupling are represented by an interac- 
tion function T which, in general, is a 27r-periodic func- 
tion of the phase difference. For weakly coupled, weakly 
nonlinear oscillators T has the universal form 0, 0] 



r(0) = (sin(0)+ 7 [l-cos(^)]) 



(2) 



Here, 0i represents the phase of oscillator i, which is cou- 
pled with strength e to a set of nearest neighbors Ni in a 



This function may be regarded as the first terms in 
a Fourier-expansion of T((jj) with the constraint that 
r(0) = 0. The symmetry breaking parameter 7 describes 
the nonisochronicity of the oscillations . 

It is well known that for sufficiently small disorder 
the oscillators eventually become entrained to a com- 
mon locking frequency fl Q, 0, Q ■ Since the oscillators 
are nonidentical, even in this synchronized state they are 
usually separated by fixed phase differences. These can 
sum up over the whole lattice to produce spatio-temporal 
patterns, which are characterized by a stationary phase 
profile, 0i — 8\. This is demonstrated in Fig. ^ where 
we have simulated system Q with interaction in a 
two-dimensional lattice. If the heterogeneity is small the 
oscillations across the lattice synchronize to a homoge- 
nous phase profile (Fig. ^1) . However, by increasing the 
disorder we observe the formation of target waves with 
decreasing wave length (Fig. QJ,-c). The emergence of 
the concentric waves is due to the symmetry breaking in 
the interaction function. For example, by reducing 7 in 
Eq. J2J the pattern becomes more irregular (Fig. ^i) . 

It is the counterintuitive observation that the isotropic 
medium with random frequency distributions of no spa- 
tial correlation (see Fig. 1 insets) can generate and sus- 
tain very regular wave patterns. Since the equations Q 
generally approximate the phase dynamics for coupled 
limit cycle oscillators this effect is not restricted to phase 
equations (JIJ . We have observed similar disorder induced 
target patterns in lattices of a variety of oscillator types, 
including predator-prey and neural Sj^stems, chemical re- 
actions and even chaotic oscillators 
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FIG. 1: Simulation results in a 2-dimensional lattice of 
150x150 phase oscillators 11121 with nearest neighbor cou- 
pling, e = 1, periodic boundary conditions, homogeneous ini- 
tial conditions and (a) 7 = 2, a = 0.029, (b) 7 = 2, a = 0.173, 
(c) 7 = 2, a = 0.433 and (d) 7 = 0.3, a = 0.433. The random 
frequencies are taken from a uniform distribution of variance 
a 2 . Plotted is the sinus of the phases #; as grey level. Similar 
results are obtained for open boundaries, (a)-(c): effect of 
increased heterogeneity; (c)-(d): influence of nonisochronicity 
7. Insets show the natural frequencies LUi as grey levels. 



In the synchronized state all oscillators rotate with 
the constant locking frequency 9i = CI, so that sys- 
tem {Q becomes a set of N equations, which have to 
be solved self consistently for the phases 6i and Cl un- 
der some imposed boundary conditions. To determine CI 
suppose first that the coupling function V is fully anti- 
symmetric r(— fa) — — r (</>), e.g. 7 = in Eq.©. In 
this case, by summing up all equations in we obtain 
Cl = uj = in the rotating frame. Thus, nontrivial lock- 
ing frequencies Cl 7^ only arise if T((f>) has a symmetric 
partr s (^) = |TO)+r(-0)) 1 



i). 



(3) 



For any coupling function T given a realization of the 
natural frequencies loi we ask for the resulting phase pro- 
file 9i. Note, that the inverse problem is easy to solve: 
for any regular phase profile 9i we can calculate Cl from 
Eq. J2J), which after inserting into Eq.QJ yields the fre- 
quencies 0J{. 

Insights into the pattern formation can be gained from 
a one-dimensional chain of phase oscillators [lrl Hij l 

n = wi + [r(&) + r(-&_i)]. 



Here we use fa = 9 i+ i — 9{ for the phase differences be- 
tween neighboring oscillators. We assume open boundary 
conditions fa — <f>N — 0. The self consistency prob- 
lem is trivial for an antisymmetric T where Cl = 0. 



In this case the T(fa) simply describe a random walk 
T(fa) — — X)}=i a; j- Thus, for small \ fa\ the phase pro- 
file 9i essentially is given by a double summation, i.e. a 
smoothening, over the disorder Ui (see Fig|2K,b). Note, 
that synchronization can only be achieved as long as the 
random walk stays within the range of T. Thus, with 
increasing system size N synchronization becomes more 
and more unlikely. 

The emergence of target waves is connected to a break- 
ing of the coupling symmetry. To explore this we study 
a unidirectional coupling with respect to <f> 



T{fa) = f{fa)Q{fa), for \fa «1, 



(•5) 



with the Heaviside function Q(fa) and f(fa>0) > 0. Here, 
the phase of oscillator i is only influenced from neigh- 
boring oscillators which are ahead of i. If the solution 
are small phase differences we are not concerned about 
the periodicity of F(</>) as the coupling is only required 
to be unidirectional close to zero. For open boundaries 
t>N = the solution to (1 lTill is given by 



»o = 



r 1 ^ 



i < m 
i > m. 



(6) 



Here, the index m is the location of the oscillator with 
the largest natural frequency, which also sets the syn- 
chronization frequency 



CI 



maxi(wi). 



(7) 



The phase differences (jSJ are positive to the left of the 
fastest oscillator, i < m, and negative to the right i > m. 
As a consequence, the phase profile has a tent shape with 
a mean slope that is given by averaging © with respect 
to the frequency distribution (FigE|c,d). We call this 
solution type a quasi-regular concentric wave. This ex- 
ample illustrates, that the asymmetry of the coupling 
function increases the influence of faster oscillators and 
effectively creates pacemakers with the potential to en- 
train the whole system. Note, that the solution Ij6l7|l is 
not possible without disorder, i.e. for a = 0. Further, 
in contrast to the antisymmetric coupling, here synchro- 
nization can be achieved for chains of arbitrary length. 

In general, the coupling function V will interpolate be- 
tween the two extremes of fully antisymmetry and unidi- 
rectional coupling in the vicinity of zero, e.g. Eq. 10 with 
7^0. As shown in Fig.2e this also gives rise to quasi 
regular concentric waves, very similar to the exactly solv- 
able system Fig. 2c. To further investigate the origin of 
these patterns note that for any given fl system Q im- 
plicitly defines two transfer maps, To : {fa_i,u>i} 1— > fa 
(4) and T^ 1 : {fa,uji} 1— > fa-i, which describe the evolution 



of the phase differences into the right or the left direction 
of the chain, respectively. The random frequencies can 
be seen as noise acting on the map (see Fig.2f) 



To 



[fi-w<-r(-&_i)]. (8) 
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FIG. 2: Phase differences (j>i (left) and phase profile 9i — 6\ 
(b,d) in units of 2n for a chain of 500 phase oscillators (QJ 
with uniformly distributed frequencies uji £ [0, 0.2] and open 
boundaries. (a,b) antisymmetric coupling @ with 7 = 0. 
(c,d) unidirectional coupling ® with r(^) = (f>Q(<f>). (e) 
System © with 7 = 2. (f) Transfer map Tn{(j>i) © for Eqs. 
@ with 7 = 2 (solid line) and fixed points <f>* (filled circles). 



It is easy to see that the breaking of symmetry leads to 
a pair of fixed points, (j)* and —</>*, in the noisy maps 



P =T n (0*,o7)=rs 1 (f). 



(9) 



geneity of the system does not permit analytical traveling 
wave solutions 6i(t) = Clt — k\i — m\. 

The general solution 1)10(1 allows for very different 
phase profiles (see Fig-EJ • The regularity of the wave pat- 
tern depends on the relative influence of the mean slope 
4>* compared to the fluctuations C,i and can be measured 
by the quality factor = <^* 2 /var(|<^ |) and the auto- 
correlation r of the (j>i. As demonstrated in Figl^both 
Qcf, and r only depend on the product 717 (see below). 
For j a -»0we find — > and r — > 1, and the solution 
is essentially a random walk (see Fig. 2a, b). With in- 
creasing values of 7(7 the correlations r are reduced and 
eventually become negative. Furthermore increases 
with the product 7CT, and for 7 > 1 can rise drastically 
(see|3t). Thus, with increasing disorder of the system 
we obtain more regular patterns until synchronization is 
lost. 

A straightforward integration of system can be 
problematic due to the long transients. During the for- 
mation of the tent-profile initially several locally synchro- 
nized clusters appear, with frequencies Qi that are deter- 
mined by the u>i in the vicinity to the cluster centers. 
Upon collision these clusters compete, where for 7 > a 
higher- frequency cluster will suppress a slower one (l2^ . 
The transient time for extinction goes with T ~ 1/Afi. 
In the disordered system with uniform p{ui) many local 
maxima with nearly identical frequencies exist, so that 
the frequency differences Ail can become arbitrary small 
with increasing system size. 

Another approach, which also applies for two dimen- 
sional lattices, relies on the Cole-Hopf transformation of 
system Assume that the fa are small so that it is 
possible to approximate the coupling function @ around 



The transfer map can be linearized at the fixed points so zer0 T(cj)) = — {e 1 ^ — l) +0((f> 3 ) After the Cole-Hopf 

transformation Oi = -h\ qi the synchronized lattice (|T|l 



and 



that Tn(0* + oj) ~ (j>* + a£ — bu; with a — r'O*) 
b = r /(0*) ■ While one fixed point, fa m the case © with 
7 > 0, is linearly stable (a < 1) the other fixed point 
is necessarily unstable. These stability properties are in- 
verted for Tq . Thus, when iterating to the right of the 
chain the fa are concentrated around fa and around — fa 
when iterating to the left. As a consequence, the general 
solution of the selfconsistency problem is build up 
from two branches around the two fixed points ±0*, su- 
perimposed by autocorrelated fluctuations d (see FigEJ) 



is reduced to a linear system 0, 0] 



4i = Eqt = 7077* -qi+ ^2 (li ~ *) 



(12) 



4j = ±0* + Ci. 



(10) 



After summation this leads to the quasi regular tent 
shape of the phase profile 6i. In linear approximation 
the Ci describe an AR(1) process 



where the random frequencies r\i — u>i/a are of zero mean 
and variance one and E = is some eigenvalue. System 
(|12f> is known as the tight binding model for a particle 
in a random potential on a lattice [2fJ. The eigenvector 
qmax corresponding to the largest eigenvalue E max will, 
in the re-transformed system of angles, outgrow the con- 
tribution of all other eigenvectors to the time dependent 
solution linearly in time. If the largest eigenvalue is non 
degenerate, the unique synchronized solution is 



Qi = ad 



(11) 



0i(*)-0i(*) = ±log(fP) 



(13) 



We want to stress that the fluctuations Q are an essential 
ingredient of the solution. Although the emerging con- 
centric waves seem to be regular the underlying hetero- 



Eq. (|13J) is well defined since the components of q*" ax do 
not change sign. Anderson localization theory |2(J pre- 
dicts exponentially decaying localized states with some 
localization length I, which after applying the reverse 
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FIG. 3: Characterization of wave patterns by the locking 
frequency -yQ (a), the wave length A/7 (b), the quality fac- 
tor Q$ (c) and the cross correlation r between neighboring 
phase differences (d). Numerical solutions are obtained by 
integrating system 11121 for one dimensional (circles) and two 
dimensional lattices (squares) in the synchronization regime. 
In one dimension the integrations were carried out with chains 
of length 500 and averaged over 50 simulation runs for 7 = 2 
(red circles) and 7 = 6 (blue circles). In the two dimensional 
system each point (green squares) represents one single sim- 
ulation in a 100x100 lattice with 7 = 2. The results using 
the eigenvector method 1121 are shown as (black +). Each 
point represents an average of 500 simulations with N = 256. 
Further, indicated in (a) (b) are straight lines with a given 
exponent a and j3 (dashed lines). The wavelength was ob- 
tained for one dimension as A = 2n/(j>* and in the two di- 
mensional system from a Fourier analysis of the phase profile. 
The plateaus in the wavelength plot are finite size effects. 



Cole-Hopf transformation yields the observed tent-shape 
phase profile with wavelength A ~ 7L Concentric waves 
emerge when A becomes smaller than the system size. 

For extremal values of 7a the system (11211 has well 
defined scaling properties E max ~ (7c) " and I ~ (7c) ^ 
[2(ij |. Perturbation theory yields a = 2 for ja <§; 1 and 
a = 1 for 7(7 3> 1. For the exponent (3 we find (3 < 1 in 
the one dimensional system, while 1 < < 2 in the two 
dimensional lattice. This implies for the synchronization 
frequency f2 and the wavelength A 



r 



A 



7 



l-Pn--f* 



(14) 



Here, 7 does not influence the wavelength as much as a 
but while an increase of 7 in one dimension leads also to 
an increasing wavelength the effect in two dimensions is 
the opposite. 

In one dimension the scaling with 7(7 holds as long as 
Eq. I|12l) can approximate the transfer map reasonably. 
The approximation breaks down when the correlation r 
of the phase differences becomes negative. In this regime 
the quality factor Q$ strongly depends on both 7(7 and 
7. The noise term bu>i in can become very small 



with increasing 7. This regime, which is not described by 
the Anderson approximation, can produce very regular 
concentric waves near the border of desynchronization. 

The constructive role of noise has often been stud- 
ied It has been shown that in spatially extended 
excitable systems noise can enhance the pattern forma- 
tion and for example is able to promote traveling waves 
[l4lll5| . Further, it is well known that disordered systems 
can synchronize faster ^{|. Here, we have investigated 
the influence of quenched noise on the pattern formation 
in the oscillatory regime. Whereas local coupling tends to 
synchronize the oscillators, the imposed disorder tends to 
desynchronize the array. We have shown that the tension 
between these two opposing forces can lead to concentric 
target patterns. 

This work was supported by the German Volkswagen 
Stiftung and SFB 555. 
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